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Abstract 

The weak coupling many-polaron formalism is applied to the case of the polaronic system con- 
sisting of impurities in a Bose-Einstein condensate. This allows to investigate the groundstate 
properties and the response of the system to Bragg spectroscopy. This theory is then applied to 
the system of spin-polarized fermionic lithium-6 impurities in a sodium condensate. The Bragg 
spectrum reveals a peak which corresponds to the emission of Bogoliubov excitations. Both ground 
state properties and the response spectrum show that the polaronic effect vanishes at large densi- 
ties. We also look at two possibilities to define the polaronic effective mass and observe that this 
results in a different quantitative behavior if multiple impurities are involved. 
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I. INTRODUCTION 



Quantum gases have revealed themselves as excellent quantum simulators for many-body 
theories from condensed matter physics, and in particular to experimentally examine strong 
coupling regimes that are not attainable with solid state experiments l| . Recently it was 
shown that when the Bogoliubov approximation is valid the system of impurities in a Bose- 
Einstein condensate Q3EC) can be added to this list through a mapping onto the Frohlich 
polaron Hamiltonianp, [J. The Frohlich solid state polaron consists of a charge carrier [elec- 
tron, hole) interacting with the LO-phonons in an ionic crystal or a polar semiconductor 0,0]. 
The Frohlich polaron Hamiltonian has resisted an exact analytical diagonalization and has 
been submitted to many approximation methods (a review on the Frohlich solid state po- 
laron can be found in Ref. j^]). The polaronic effect of an impurity in a BEC has been 
the subject of several recent theoretical studies. Some examples are the polaronic effects in 
optical lattices^- 9], the application of the Feynman variational path integral technique to 
study the ground state properties [101]. an extension of this technique to examine the response 



properties 



111 ] and a strong coupling approximation 12j. This revealed that the static prop 



erties are characterized by the ratio of the masses of the bosons and the impurity together 



with the polaronic coupling parameter a, defined as: 



a 



"'IB 

£,o>bb 
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with a/sthe impurity-boson scattering length, <ibb the boson-boson scattering length and 
£ the healing length of the condensate. As a function of this cou plin g parameter two pola- 
ronic regimes were identified, reminescent of the acoustic polaron [13(. In the context of the 
Frohlich solid state polaron the intermediate and strong coupling regimes still supply im- 
portant theoretical challenges but the weak coupling regime is well understood. The weak 
and intermediate coupling theory of Lee, Low and Pinesjl^] was extended by Lemmens, 



Brosens and Devreese 



15] to the case of an interacting polaron gas. Similarly, the optical 



n 



absorption theory of Devreese, Goovaerts and De Sitter [l6| was extended to the case of in- 
teracting polarons by two of the present authors, in Ref. 17J . In these extensions, the effect 
of the interactions between the polarons is taken into account through the structure factor 
of the electron gas. The goal of this paper is to apply a similar generalization to extend 
the polaronic theory of a single impurity in a Bose condensate to the case of a dilute gas of 
interacting impurities. This is needed since an experimental realization of the BEC-impurity 



polaron will involve multiple impurties and will require the use of Feshbach resonances to 
reach the strong coupling regime (iq|. 

Within the Bogoliubov approximation the Hamiltonian of N impurities in a condensate 
can be mapped onto the sum of the mean field energy Emf and the Frohlich iV-polaron 



Hamiltonian H^ ol 



10j. The mean field energy is given by: 



E MF = E GP + N NU IB [q = 



(2) 



where the first term is the Gross-Pitaevskii energy of the condensate 18[] and the second 
term is the interaction shift due to the impurities with Nq the number of condensed bosons 
and Uib {q) the Fourier transform of the impurity-boson interaction potential. The polaron 
Hamiltonian describes the mutual interaction between the impurities and their interaction 
with the Bogoliubov excitations: 



N 



N 



k k 1 



r | /■, — Tj \ . (3) 



The first term represents the kinetic energy of the impurities with mass mi and position 
(momentum) operators (pj, the second term gives the kinetic energy of the Bogoliubov 
excitations with dispersion uj^ and creation (annihilation) operator at (%), the third term 
is the interaction between the impurities and the Bogoliubov excitations with interaction 
amplitude Vt and the fourth term represents the mutual interaction between the impurities 
with v (f) the interaction potential. The Bogoliubov dispersion is given by: 



(4) 



where the speed of sound in the condensate was introduced: c — h/ (v^mB^) , with tub the 
mass of the bosons. The interaction amplitude is given by: 

1 1/4 

(5) 



Vt 



'N 



9 IB} 



.(#) +2. 

where a contact potential was assumed for the impurity-boson interaction potential: 
Uib (q) — Qib- At low temperature the amplitude gjB is completely determined by 
the effective mass m r = (m^ 1 + mj 1 ) 1 and the impurity-boson scattering length clib' 
giB = 2nh 2 a IB /m r . 
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Another important consideration is the stability of the mixture against phase separation. 
At zero temperature it was shown in Refs. [3] and [^(J that the mixture will be stable if 
the following inequality is satisfied: 

1/3 < (6vr) 2/3 m r a BB 



127T mBfriia 2 IB ' 

with n the density of the impurities. 

Since the impurities are not charged the response can not be studied through optical 
absorption measurements as for the Frohlich solid state polaron. In Ref. ll| it was shown 
that the internal excitation structure of the BEC-impurity polaron can be probed with 
Bragg spectroscopy which is a technique that has proven to be very successful in the study 
of Bose-Einstein condensates (see for example Refs. 21] and 22]). The experimental setup 
consists of two laser beams with wave vectors k\ and hi and energies U\ and U2 which are 
impinged on the impurities. These can then absorb a photon from one beam and emit it in 
the other beam which results in the exchange of a wave vector h = h\ — hi and an energy 
fkd = hu\ — hu2 to the impurities. The response of the system after an exposure during a 
time interval r can be measured by counting the number of atoms Nsragg that have gained 
a wave vector k as a function of u> which in the formalism of linear response theory is given 
byQ: 

2 fV x 2 



N Bragg = JI J T Im X (k, U J , (6) 

with V the amplitude of the laser induced potential and x (k> w J the density response 
function, defined as: 

In the following we start by summarizing the main results of Refs. [ljj] and [l^] on 
the weak coupling treatment of the Frohlich solid state iV-polaron system are summarized 
together with an indication of what changes if one considers impurities in a BEC. We also 
show two different ways to define the polaronic effective mass. This formalism is then 
applied for spin-polarized fermionic impurities and the results are examined as a function 
of the impurity density and the exchanged momentum for the Bragg response for lithium-6 
impurities in a sodium condensate. 
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II. WEAK COUPLING TREATMENT OF THE MANY-POLARON GAS 



A. Ground state properties of the many-polaron gas and the inertial effective 



mass 



The weak coupling variational method as introduced in Ref. 14] by Lee, Low and Pines 
for the description of a single polaron was generalized in Ref. (hi]] to the case of a many- 
polaron system. With this purpose the following variational wave function was introduced: 



\^LDB >=U\0) 

with |0) the vacuum wave function for the bosonic excitations, 
impurities and U a canonical transformation of the form: 



N 



are 



ik.r 
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the wave function of the 



(9) 



where {fj:} are variational functions. Minimizing the expectation value of the iV-polaron 
Hamiltonian ((3]) with respect to the variational wave function (jSJ) as a function of {f^} 
results in the following expression for the ground state energy: 



S k 



S 2 k) V; 



i ^kS[k) + r S 



where e^ in is the kinetic energy per particle: e^i 
static structure factor of the impurities: 



N- 1 (if. 



^i 2mj 



(10) 



ip ) and S ( k ) is the 



S k 




N 



ik.(fi-rj) 



(ii) 



Another important property of the polaronic system is the effective mass of the polarons 



m* . Within the one-polaron weak coupling formalism m* was calculated in Ref. [14J. The 
generalization to many polarons allows a determination of this effective mass which is related 
to the many-polaron system as a whole and shall be called the inertial effective mass in the 
following. The total momentum V =Y2iPi + Ylk^^l^k commutes with the Hamiltonian 



dSJ) and thus is a conserved quantity. This conserved quantity can be explicitly introduced 
in the minimization process by means of a Lagrange multiplier v, which corresponds to the 
velocity of the polaronic system: 

Ki (*o = Ki - * \Tpi + E ^4% - A • ( 12 ) 

minimization of the expectation value of ( I12p with respect to ([H]) as a function of {/^} and 
v, together with an Taylor expansion for small v, results in: 

2 , \V t \ 2 S 2 (k) 

r = N mi v + -h 2 Nj2 — — — — — s* (13) 

3 s fas(k) + §£) 

where isotropy of the system was assumed. We can now identify the inertial effective mass m* 
through the expression V = Nm*v, which relates the speed of the polaron to its momentum, 
i.e.: 



m = mi 



3 s (%SW+S0 



Note that the inertial effective mass (I14p is a property of the whole polaron system. This 
is important for an experiment where the behavior of all the impurities together is studied. 
For example when one studies the collective oscillation of the impurities in a harmonic trap 
the resulting effective oscillation frequency is a function of this inertial effective mass. This 
effect has been used to experimentally determine the effective mass of Fermi polarons in 
Ref. |23]. 

Note that for the ground state properties the presence of multiple impurities is completely 
described through the static structure factor of the impurities. In the limit of vanishing 
density the static structure factor becomes 1 and the expressions for the energy and the 



effective mass of a single polaron from Ref. 



14( are retrieved. 



B. Bragg response of the many-polaron gas and the spectral effective mass 



In Ref. 17] the optical absorption of the iV-polaron gas consisting of electrons interacting 
with phonons was calculated within the weak coupling formalism. The polaronic system 
consisting of impurities in a BEC can not be probed with optical absorption but instead an 



appropriate experimental technique is Bragg spectroscopy as shown in Ref. 



ll| . This means 



that a finite momentum exchange has to be taken into account which is negligible in the 
case of the optical absorption calculation. Furthermore it is the density- density correlation 
function that determines the response instead of the current-current correlation which is 
needed for the optical absorption. These two are however closely related since the application 
of two partial integrations transforms equation (j7J) into a correlation function of the time 
derivative of the density after which one can use the Fourier transform of the continuity 
equation (pg = ik.j£) to obtain a current- current correlation function: 



X [u,k 



hw 2 



dte 



iujt 



k-j% (t) , k.j 



For the calculation of (1151) we use the derivation of Ref. 



(15) 



17] which is based on the wave 
function (jSJ) and which preserves only terms of lowest order in the coupling amplitude \VA . 
The main difference of the present result is the incorporation of a finite momentum exchange 
k. This results in the following expression for the imaginary part of the density response 
function (TT5T) which is proportional to the Bragg response ([6]): 

lm X = ^j^2 (^) \V/s(q + k,u-u^ , (16) 

1 q 

where S (k, uij is the dynamic structure factor of the impurities: 

1 



S [k, u 



2nh 
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(17) 



The response also allows a determination of an effective mass which we shall call the spec- 
tral effective mass. This is done through an extension of the f-sum rule that was introduced 
in Ref. 24] and generalized for the polaron in a condensate in Ref. 11] to: 



h / au uj lim 

2m* J + k^ 



Im 



X (k,u 



Ntt 



P 2m/ (18) 

The origin of this sum rule is a £-peak at u = in the response, the spectral weight of this 

5-peak in the k — > limit equals the first term in (j!8[) . It is important to note that the 

spectral effective mass as determined by (I18p is a single particle property and is thus not the 

same as the inertial effective mass of expression (|14p . Only in the limit of one polaron the 

two masses coincide as can be checked from equation (1181) . This polaronic spectral effective 

mass, as defined in (|18|) . was for example measured in Ref. [25] for the Frohlich solid state 

polaron. 

Note that in this formalism the influence of multiple impurities on the polaronic response 
is solely determined by the dynamic structure factor of the impurities. 



III. SPIN-POLARIZED FERMIONIC IMPURITY GAS IN A BEC 



In this section we apply the results of the previous section to the case of a spin-polarized 
gas of fermionic impurities. Because of the Pauli exclusion principle the s-wave scattering 
length is zero which results in essentially no interactions between the impurities at low 
temperature. This means the impurities can be described as an ideal fermionic gas for 
which the static and dynamic structure factor are known. All results in this section are 
presented in polaronic units, i.e. h = mi = £ = 1, and for lithium-6 impurities in a sodium 
condensate, i.e. m B /mj « 3.82207. 



A. Ground state properties 



The static structure factor of an ideal fermionic gas at zero temperature is given by 26] : 



3 ^ - ~ (ttM 3 if k < 2k f 



S [k) = { 22k - *\»») F , (19) 

1 if k > 2k F 

where kp is the Fermi wave vector which, for a non-degenerate gas, is given by kp = 

1 /3 

(6ir 2 n) , with n the impurity density. Introducing the kinetic energy e^m of an ideal 
fermionic gas together with the Bogoliubov dispersion (jlj) and the interaction amplitude ([5]) 
in the expression for the ground state energy (Tl0]) leads to: 



mBS 2 [k) / k 2 \i/a 



m B 



AT 10 2tt V m B J J \ji^W+2S{k)+m B k 2 \ k2 + 2 J rn B + l 

(20) 

with a the polaronic coupling parameter from expression ([1]). The first term represents 
the kinetic energy and the second is the polaronic contribution, the interaction energy has 
vanished since we are describing an ideal gas. Notice that in the polaronic contribution an 
additional term has appeared which is needed to obtain a convergent energy and was also 



obtained in Ref. 10j as a result of the renormalization of the boson-impurity interaction. 
The polaronic contribution in equation (|20|) grows linearly with kp as the number of particles 
is increased. Since the kinetic energy of the Fermi gas grows as k F , the relative contribution 
of the polaronic energy with respect to the kinetic energy decreases. In the limit of large 
densities, the kinetic energy dominates the polaronic effect. 
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FIG. 1: The inertial effective mass of the polarons as a function of the density for lithium-6 
impurities in a sodium condensate at a = 0.01 according to formula (|2ip . The inset shows the 
behavior at small densities. In the limit n — > the effective mass of a single polaron is retrieved 
and for n — > oo the effective mass becomes the bare impurity mass. 




Introducing the Bogoliubov dispersion (j3j) and the interaction amplitude (jSJ) in the ex- 
pression for the polaron inertial effective mass (I14p results in: 

4a (m B + l\ 2 r „ 1 k 2 S 2 (k 
m} = 1 + — / rffc ^ 3. (21) 
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This expression is presented in figure [T] as a function of the impurity density. In the limit 
n — > the one-polaron result is retrieved, which was already anticipated using formula (|T4|) . 
For n — > 00 the inertial effective mass becomes equal to the bare impurity mass. This limit 
can easily be examined analytically with equation (12 ip which reveals a n _1//3 behavior at 
large densities. For the Frohlich solid state polaron the same qualitative behavior as in figure 



[T]was found for the effective mass if the electrons are described as a free gas 27|, 128 ]. 



B. Response to Bragg spectroscopy 



For the response the dynamic structure factor is needed, which for an ideal fermionic gas 
at temperature zero is given by |29l |: 

2 



S (q, u) 
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UJ 



(22) 
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FIG. 2: The Bragg response (123ft of the polaronic system consisting of polarized lithium-6 impurities 
in a sodium condensate as a function of to for different values of the exchanged momentum k. The 
impurity density is taken n = 0.01 and the polaronic coupling parameter is a = 0.01. 

with 9 the Heaviside step function. Introducing the Bogoliubov dispersion ((4]) and the 
interaction amplitude ([5]) in the expression for the imaginary part of the density response 
function ( IT6|) leads to: 

In figure [2] expression (1251) is shown as a function of u for various momentum exchanges. 
A peak is seen that represents the emission of Bogoliubov excitations and which is shifted 
to higher frequencies for larger momentum exchange. This behavior is to be expected since 
more energy is needed to create a Bogoliubov excitation with a higher momentum. 

In figure [3] the peak in the Bragg response is presented for different impurity densities. 
Note that for higher densities the spectral weight of the u > peak diminishes. This is not 
in contrast with the f-sum rule: the spectral weight of the u = delta peak in the spectrum 
compensates, in accordance with expression ffTSj) . For k — > the weight of this 5-function 
is related to the spectral effective mass such that an attenuation of the peak at u > 
corresponds to a decrease of the spectral effective mass. The spectral effective mass deduced 
with the sum rule (|T8|) is presented in figure H] as a function of the impurity density. Note 
that the spectral effective mass behaves generally only qualitatively the same as the inertial 
effective mass from figure [T] and only in the limit n — > the same result is retrieved. This 
is because the spectral effective mass is a one-particle property while the inertial effective 
mass is related to the entire polaronic system. 

We would like to emphasize that according to our calculations different definitions of the 

10 
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FIG. 3: The Bragg response (I23j) of the polaronic system consisting of polarized lithium-6 impu- 
rities in a sodium condensate as a function of uj for different impurity densities n. The exchanged 
momentum is taken k = 1 and the polaronic coupling parameter is a = 0.01. 
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FIG. 4: The spectral effective mass of the polarons as a function of the impurity density for 
lithium-6 impurities in a sodium condensate at a = 0.01 as determined by the sum rule (|18p . 



effective mass result in another behavior and can in general not be compared. This also 
means that for an experiment it is important to know which effective mass is of importance 
for the specific setup. 



IV. CONCLUSIONS 



The weak coupling many-polaron formalism which was developed in the context of 



Frohlich solid state polarons in Refs. 



151 ] and [17| was applied in the present work to 



the case of the polaronic system consisting of impurities in a Bose-Einstein condensate. The 
properties of the ground state and the response to Bragg spectroscopy were examined. The 
sum rule of Ref. 24j which relates the response to the spectral effective mass was formulated 
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in the present context. Also a calculation of the inertial effective mass was presented which 
is related to the many-polaron system as a whole. It turns out that the ground state prop- 
erties are determined by the static structure factor of the impurities while the response is 
governed by the dynamic structure factor, as in the case of the Frohlich solid state polarons. 

This generalization of the many-polaron formalism was then applied to the case of spin- 
polarized fermionic impurities which behaves as an ideal gas for which the structure factors 
are well-known. The numerical calculations were done for lithium-6 impurities in a sodium 
condensate. Both the ground state properties and the Bragg response indicate that in 
the limit of high impurity density the polaron effect disappears, which is also the case for 
Frohlich solid state polarons. In the Bragg response a peak was observed that corresponds 
to the emission of Bogoliubov excitations. The behavior of this peak as a function of the 
exchanged momentum and the impurity density was examined. It was also shown that the 
two definitions for the effective mass exhibit a different behavior, this is because the inertial 
effective mass is a property of the system as a whole while the spectral effective mass is a 
one-particle property. 
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